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We reconsider spherically symmetric black hole solutions in Einstein-Aether theory with the con- 
dition that this theory has identical PPN parameters as those for general relativity, which is the 
main difference from the previous research. In contrast with previous study, we allow superluminal 
propagation of a spin-0 Aether-gravity wave mode. As a result, we obtain black holes having a 
spin-0 "horizon" inside an event horizon. We allow a singularity at a spin-0 "horizon" since it is 
concealed by the event horizon. If we allow such a configuration, the kinetic term of the Aether 
field can be large enough for black holes to be significantly different from Schwarzschild black holes 
with respect to ADM mass, innermost stable circular orbit. Hawking temperature, and so on. We 
also discuss whether or not the above features can be seen in more generic vector-tensor theories. 
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I. INTRODUCTION 

Identifying the contents of dark energy and dark mat- 
ter (DE/DM) is one of the most important subjects in 
cosmology. It is frequently argued that gravitational the- 
ories are an alternative to DE/DM. Recently, tensor- 
vector-scalar (TeVeS) theories have attracted much at- 
tention since they do not only explain galaxy rotation 
curves but also satisfy many constraints from solar ex- 
periments Although a deficency in explaining the 
mismatch between luminous and dynamical masses in 
clusters of galaxies by TeVeS has been pointed out 0, 
resolution of this problem by considering a generalized 
vector-tensor theory has also been reported [3]. More- 
over, these vector fields might explain an accelerated ex- 
pansion of the universe [3Uj| and might be important in 
inflationary scenarios [aly]- The origin of such a vector 
field is argued in ff\. 

However, it is nontrivial whether or not these theories 
satisfy the constraints by strong gravity tests. Notice the 
result for scalar-tensor theories where compact objects 
have strong deviations from those in general relativit 
(GR) even in the cases that satisfy weak field tests 
To study vector fields in a general form is difhcult. For 
example, results in TeVeS are still limited to special cases 
such as 0. Thus, as a first step, it is important to inves- 
tigate a simplified model which is tractable and instruc- 
tive for general cases. One such useful model would be 
Einstein- Aether (EA) theory [l^], where all parameter- 
ized post-Newtonian (PPN) parameters [llj can be the 
same as those in GR EA theory is a vector-tensor 
theory, and TeVeS can be written as a vector-tensor the- 
ory which is the extension of EA theory [l^. In EA 
theory, strong gravitational cases includin g b lack holes 
have been analyzed to some extent [H, flsLIl^ [TtI [Tsj . 

Nevertheless, the analysis of black holes has been lim- 
ited to the case in which the event horizon coincides with 



the spin-0 horizon [15], and this case does not neces- 
sarily satisfy weak fields tests. Thus, it is interesting 
to ask whether or not significant differences from the 
Schwarzschild black hole appear when weak fields tests 
are satisfied. For this reason, we argue black holes with 
the case in which the EA theory has identical PPN pa- 
rameters as in GR. 

This paper is organized as follows. In Sec. II, we ex- 
plain EA theory and summarize constraints located by 
previous research. In Sec. Ill, we mention our method 
of analyzing black holes. In Sec. IV, we show the results 
and compare them with the Schwarzschild black hole. In 
Sec. V, consequences and future subjects are discussed. 
In Appendix A, we summarize basic equations. We use 
units in which c — 1 and follow the sign conventions of 



Misner, Thorne, and Wheeler 
metrics. 



e.g., (-,+,+,+) for 



II. EINSTEIN-AETHER THEORY 



A. The action and basic equations 



We consider the following action [Tl 
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where u° is a vector field and := u°'Ua- Ci 
1,2,3,4) are theoretical parameters in EA theory. A is 
a Lagrange multiplier ensuring the vector field to be 
unit timelike vector everywhere. 

Varying this action with respect to A and we have 
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Multiplying Eq. (12. 5p by Ua, we have 

A = C4U^ + u''V„J" . (2.8) 
Varying the action with respect to the metric, we have 



Gab = Vr 



where 



J (a^b) - J(a "b) + ^ (ab)1 



+Cl (VaUmVfcU™ - V„WaV"wb) 
+C4?iaUfc + \UaUb - -Qab^n , (2.9) 
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where we eliminate C2 and C4 with Eq. (|2.12p . 

For these sound velocities to be equal to or larger than 
the photon velocity, or, to ensure stability against linear 
perturbation in Minkowski (or FRW) background and 
linearized energy positivity, we have 0, [H, [2J, |2^ 



< c+ < 1 , 
< c_ := ci 



C3 < 



3(1 -c+) 



(2.18) 



Radiation damping was also analyzed in [26|, |27| , which 
almost restricts function of c_ based on the ob- 

servation of, say, B1913-I-16 ^28i] . 



B. Present Constraints in EA Theory 



III. ANALYSIS IN A SINGLE-NULL 
COORDINATE SYSTEM 



If we assume the weak field and slow-motion limits in 
EA theory [l3|, we have to take Newton's gravitational 
constant as 



Gn 
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(2.11) 



to reproduce Newtonian gravity correctly. For all the 
PPN parameters to coincide with those in GR, we have 



C2 



-2ci - C1C3 -I- C3 
3ci 



Ci 



(2.12) 



If we assume Friedmann-Robertson- Walker (FRW) 
space-time and the Aether is aligned with a cosmolog- 
ical rest frame, the cosmological gravitational constant is 
given by [13] 
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where c+ :— ci + C3. Using primordial *He abundance, 
we have 



|G, 
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II < 1/ 



(2.14) 



From the maximum mass of neutron stars ^ 2Mq |21| . 
[2^ . we have ci -I- C4 < 0.5 ^ 1.6, depending on EOS [l8l |. 

In [2^, the sound modes are analyzed by expanding 
the metric and the Aether around the Minkowski met- 
ric. As in the case in GR, we have two spin-2 modes. 
As peculiar to EA theory, there are three wave modes. 
Two correspond to a transverse spin-1 mode, and one 
corresponds to a longitudinal spin-0 mode. The squared 
speeds of them are summarized as 



Cl3 



3(ci - C3)(l - C13) 

ci(2ci - cl + cl) 

2(1 - Ci3)ci3(ci - C3) 
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(2.15) 
(2.16) 
(2.17) 



Our purpose in investigating black holes in EA the- 
ory is not to give a further restriction but to understand 
generic features of vector-tensor theories under the con- 
dition that weak gravity tests are satisfied. This is the 
main difference from the previous research Il5|, which 
investigates black holes with the parameters [29j 



C2 



-(Cl + C3)^(ci - C4) - 2(C3 + C4) 

(ci — C4)(3c3 -I- 4c4 — Cl) + 2 



(3.1) 



and C3 — — C4, or C3 — — 2c4 -I- ci, or C3 = — ci. In these 
parameters, qualitative differences from Schwarzschild 
black holes have been shown. It is nontrivial whether 
or not this is true even for the case which satisfies weak 
gravity tests. 

From this point of view, we take the following strategy, 
(i) We assume (j2.12|) since the constraints by the solar ex- 
periments are severe, (ii) We assume (I2.18p . Otherwise, 
a naked singularity appears outside the event horizon in 
general. As for other constraints, notice that (|2.14p is 
satisfied if (|2.12p is satisfied. Constraints from neutron 
stars and from radiation damping are related to strong 
gravity tests at least partially. For the above reasons, 
we do not impose these constraints. Thus, we have two 
theoretical parameters (c+, c_) with the condition (|2.18p . 

We write a static and spherically symmetric line ele- 
ment in a single null coordinate system as. 



'N{r)dv^ + 2B{r)dvdr + r^dfl^ 



(3.2) 



In this coordinate, the vector field takes the form of 

u = a{r)dy + b{r)dr ■ (3.3) 

b{r) 7^ means that the Aether is not aligned with the 
timelikc Killing field, which is inevitable because of the 
event horizon. From Eq. 



Na^ + 2Bab 



-1 . 



(3.4) 



We can eliminate A with Eq. (|2.8p . Then, from the Ein- 
stein and Aether equations, we obtain basic equations. 
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which can be written schematically as 

N' - h{B,N,a,a') , (3.5) 
B' = f2iB,N,a,a') , (3.6) 
a" - MB,N,a,a') , (3.7) 

where the prime denotes the derivative with respect to r. 
Here, we have eliminated b with Eq. (|3.4p . The explicit 
form is summarized in Appendix A. 

The boundary condition at the horizon rh is N(r\^) — 0. 
We set B{r\y) = 1. We can also set rh = 1 since there is 
no scale in the present theory. In this sense, it is assumed 
that the area coordinate r is normalized by the horizon 
radius below. 

If we use a rescaling freedom of v as dv' — B{oo)dv, 
the asymptotic form of the metric is written as 

iV(oo) 



(a) 



rdv'^ + 2dv'dr + r^dn^ 



(3.8) 



B(00)2 

Thus, the boundary condition at spatial infinity for the 
asymptotic flatness is 

A^(oo) = B{oof . (3.9) 

We should require 

&(oo) = , (3.10) 

for the Aether to be aligned with the timelikc Killing 
field. Then, by Eq. ([^ . we have 

1 



a{oo) 



(3.11) 



B{oo) 

We can determine the pair of ah := a(?'h) a-nd aj^ :— a'(rii) 
as shooting parameters, one of which is fixed by (|3.1ip . 
Thus, there remains one freedom. Fixing this freedom is 
done as follows. 

Even in the spherically symmetric case, there is a spin- 
mode. Then, we can define the effective metric for a 
spin-0 mode as 



5ab-[(so) ~l]UaUb. (3.12) 

We call the horizon associated with this metric as the 
spin-0 horizon. The freedom mentioned above is fixed by 
the requirement that the regularity at the spin-0 horizon 
which is inside the event horizon. 

However, since the asymptotic observer is insensitive 
to the regularity at the spin-0 horizon, we permit the sin- 
gularity at the spin-0 horizon. For this reason, we leave 
one freedom. In concrete terms, we obtain Oh iteratively 
for some a'-^, which is regarded as a free parameter. We 
use the Bulirsch-Stoer method in our numerics l30l. 



IV. PROPERTIES OF SOLUTIONS 



Mass and Hawking temperature of EA black 
hole 




OQ 1 



(b) 




We show several asymptotically flat solutions in Figs.[T] 
(a)-(c) for c+ = 0.4 and c_ — 0.1. In the figures, we have 



FIG. 1: Field configurations for c+ = 0.4 and c_ — 0.1. De- 
noted numbers in each figure, ranging from —1 to 1, represent 
the values of aj^. We normalize the quantities Gm and r by 
the horizon radius r^. The solution with the smallest a{j has 
largest deviation from a Schwarzschild black hole. 
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we define the mass function r7i(r) by 
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FIG. 2: (a) a!^ v.s. GMadm and (b) v.s. GMtot for 
several sets of c+ and c_ . Physical quantities are normalized 
by the horizon radius rh. Notice that there is a lower limit 
'Jh.crit below which there is no regular solution. Near a(j ^.^it, 
GMadm and GMtot depend on a(j remarkably. 



selected five solutions. The differences of these solutions 
are the changing boundary value af^, ranging from —1 to 
1, as denoted in the figures. Figure [1] (a) shows that we 
can determine an ah that satisfies the asymptotic condi- 
tion (|3.1ip for various values of a{j. We also show B{r) in 
Fig. [T] (b). Since B{r) = const. = 1 for a Schwarzschild 
black hole, it indicates that there are differences in phys- 
ical quantities from those for Schwarzschild black holes. 
Figure [T] (c) shows a "mass" function. In AE theory, it 
is important to distinguish different notions of mass. If 



m(r) 



r 

2G 



1 - 



N 
52 



(4.1) 



we can interpret m(oo) as ADM mass Madm- As we 
can see, m(r) monotonically decreases. Our calculation 
suggests that this is generic. This is not surprising since 
energy conditions are not necessarily satisfied in EA the- 
ory |25j . 

Since Figs. [1] show that the deviation from the 
Schwarzschild black hole is largest for the smallest value 



of, 



'h' 



it is natural to ask whether or not there is a lower 



limit flj'j below which there is no regular solution. We 
show the relation aj^ and Mabm for various values of 
c+ and c_ in Fig. [2] (a). Typically, Madm is smaller 
than that of a Schwarzschild black hole by about 10%, 
which is consistent with the result in [lB|. Here, we ob- 



tain Oh iteratively to satisfy Eq. 



for each 



For 



a , < a' 



h,crit 



jj, we cannot find an appropriate value of a^- 



^h,crit 



depends on c_|_ and c_. As approaches a'^ 



crit ' 



dMji^DM/da'^ tends to diverge. Since we obtain solutions 
numerically, it is nontrivial whether Madm is bounded 
or not from below. In particular, it is important to re- 
veal the positivity of Madm. However, since the energy 
conditions are not guaranteed [25l | , we cannot prove it at 
present. 

For Madm, the difference caused by the change of 
c_ is not clear. We can define total energy Mtot by 
GnMtot = GMadm since the gravitational constant we 
feel is different from that in GR as seen in Eq. (|2.1ip . 
We also exhibit the relation oJj-Mtot in Fig.[2](b). This 
figure shows the differences caused by the change of c_ . 
Mtot decreases as c_ increases as similar to c+. 

If we contemplate these diagrams from a different view- 
point, we notice that the horizon radius of black holes in 
EA theory is larger than that of a Schwarzschild black 
hole for fixed GMtot (or GMadm)- Therefore, one 
might think that black holes in EA theory have larger 
entropy. However, since we have the Lorentz violating 
field, it is nontrivial to establish black hole thermody- 
namics [Sll, [12] ■ Thus, the comparison of the black hole 
entropy, which is crucial to discuss the stability of black 
holes, belongs among our future tasks. 

It is also important to reveal what happens at the 
critical point, ~ ^hcrit- "^^^ '^^y Point is the factor 
(—1 -t- Na'^) in the denominator in (|A4[) . For — aj^ 
(—1 + Na^) becomes zero at finite r. Thus, solutions 
disappear. We also show the dependence of Hawking 
temperature Th for c+ = 0.6 and c_ =0.1 in Fig. [S] 
From this diagram, it is supposed that Th diverges for 



the solution at ai 



-'h,crit ■ 



Thus, It is intriguing to 



consider an evaporation process of such black holes. 
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FIG. 3: Hawking temperature Th (normalized by rh) for c+ — 
0.6 and c_ =0.1 suggesting that Th diverges for the solution 

S-t Oh.crit- 

B. ISCO of EA black hole 

We shall turn to more realistic problems. We consider 
the possibility of distinguishing black holes in EA the- 
ory from Schwarzschild black hole by observation. In 
Ref. fl^, the innermost stable circular orbit (ISCO) for 
neutron stars in EA theory was analyzed. The result is 
that the deviation from the Schwarzschild black hole is at 
most several percent. But this is not necessarily the case 
in the present situation, as shown below. The differences 
occur since we have the freedom parameterized by and 
the Aether is not static. These facts will be important 
if we consider observations such as Constellation-X [ssj . 
which tracks the motion of individual elements near black 
holes. 

From an equation for timelike geodesies for a unit mass 
particle, we have effective potential V as 

y(r) = ^ + l) . (4-2) 

where L is the angular momentum normalized by the 
horizon radius. 

We show the typical configurations of V in Fig. 2] for 
EA theory (with c+ = 0.6, c_ = 0.1 and aj^ = 0.78 ~ 
a(j ^j.jj) and for GR (Schwarzschild black hole), where the 
angular momentum of the test particle L is fixed as i = 
1.5. We find that a potential minimum exists even for 
L = 1.5 in EA theory. 

We show the dependence of nsco (normalized by rh) 
on a(j in Fig. [5] (a). Notice that nsco = 3 for the 
Schwarzschild black hole. Therefore, the difference is 
nearly 10% for a(j ~ Ohcrif ^^^o impressive to 

write the ISCO normalized by GMtot (or GMadm), 



which is shown in Fig. [5] (b). In this case, we can 
find the difference from the Schwarzschild black hole 
(hsco/GMadm = 6) is more than 20%. 

Finally, let us comment on the parameter region of 
(c+, c_) in which black hole solutions exist. We obtained 
solutions even for c+, c_ > 1, which seems to conflict 
with the previous results [15]. However, since we do not 
assume regularity at the spin-0 horizon against the case 
in [l5[, it is not inconsistent. The qualitative properties 
are same as in other parameter regions, although quanti- 
tative differences from Schwarzschild black hole become 
larger for large (c+ , c_ ) as we expect from Fig. [51 These 
features are same as in [l5| where the consistency with 
the weak gravity tests are not necessarily imposed. 



> 




FIG. 4: The potential V for EA theory (c+ = 0.6, C- =0.1 
and a{j = 0.78 ~ a{j ^^^f.) and for a Schwarzschild black hole 
where the angular momentum of the test particle L (normal- 
ized by Th) is fixed by L = 1.5. There is a potential minimum 
in EA theory while there is none for a Schwarzschild black 
hole. 



V. CONCLUSION AND DISCUSSION 

We have reanalyzed black hole solutions in EA theory 
while assuming that all the PPN parameters are the same 
as those for GR, resulting in two theoretical parameters 
C-i- and c_. This is a main difference from the previous 
study As another difference, we do not assume reg- 
ularity at the spin-0 horizon since this is inside the event 
horizon. Interestingly, we find a(j ^^.j^ below which there 
is no regular black hole solution. Near ajjcrif 
viation of black hole mass Mtot (or Madm) and ISCO 
risco from those for the Schwarzschild black hole become 
large. 

These results are instructive for other cases. If we con- 
sider the case with rotation, freedom of the vector field 
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(a) 




(b) 




FIG. 5: The ajj dependence of the innermost stable circular 
orbit (ISCO) for EA theory with c+ = 0.6 and c_ = 0.1. (a) 
ISCO normalized by rh. (b) ISCO normalized by GMtot 
and GA/adm. 



is added to (|3.3p . Then, it also contributes the kinetic 
term of the vector field, enhancing the differences from 
the vacuum solution. This would also be true in other 
vector-tensor theories. For this reason, it is important to 
consider rotational black holes in vector-tensor theories, 
if we are to constrain them. 

Although we have revealed many properties of EA 
black holes, some important problems remain to be in- 



vestigated. One is the positivity of the energy, which is 
necessary for the stability of the system. Related to this, 
to establish the black hole thermodynamics is also im- 
portant. As a consistency check, we should also perform 
the linear perturbation for the black holes [SJ]. 

The other is whether or not regular spin-0 horizon hap- 
pens as a result of gravitational collapse. In it is 
shown that regular spin-0 horizon happens if we consider 
a gravitational collapse of a massless scalar field. Thus, it 
is important to investigate this feature in a general case. 
It is also interesting to investigate the critical behavior 
of such a system 36] . Of course, these are not problems 
particular only to EA theory but also issue confronting in 
more generic vector-tensor theories. Thus, it is desirable 
to investigate them in a unified way. 
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APPENDIX A: BASIC EQUATIONS FOR 
EINSTEIN-AETHER SYSTEM 



The equation for N is 



J2 Hia' = 



(Al) 



where 



He = c+[(3c_ - c+)7V2 + 2(3c_ -I- c+)rNN' + 
(3c_ - c+yN''^ , 

H5 ^ 2c+rA^[(3c_ + c+)N + (3c_ - c+)rN']a' , 

H4 = -12(c_ + c+)B(rf + 2[-3c_(-2 + c+) + 

c+(6 + c+)]N ~ 2[3c_(-2 + c+) + (-6 + c+)c+]rN' + 

(3c_ - c+)c+r'^N'^a'^ , 

H3 = -2c+{3c^+c+yN'a' , 

H2 = -c+a2[-3c_ + C+ + 2(3c_ + c+yNa'^] , 

Hi = -2c+(3c_ + c+)ra' , 

Ho = (3c_ - c+)c+r^a'^ . 

(jAip is the quadratic equation for N' (Notice Hq). If 
we solve (jAll) about TV', we obtain the equation which 
satisfies asymptotically flatness as 



(3c_ 



(A2) 
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where 

hs = -c+(3c_ + c+)N , /i2 = c+(c+ - 3c_)rA^a' , 

= 3c_(c+ — 2) + (c+ — 6)c+ , /lo = c+(3c_ + c+)m' , 
X = [-3c2_(c+ - 1) - (c+ - 3)4 + c_c+(6 - 4c+ + 
cl)]a'^ + c+[(3c?. + 2c_c+ - cl)B^ + 2c+(c_ + c+ - 
C-C+)N]a'^ + c-^iV^a^ + c+[3ci(c+ - 1) + 
c_(c+ — 4)c+ — c'^]raa' — c+[3c^(c+ — 1) + 
c_(c+ - 2)c+ + c+]r7Va^a' + c^c^r'^a''^ . 

Notice the denominator in (|A2|1 . For c+ — 3c_ , we should 
use (jAT|) . 

The equation for _B is 



c_4(13 + 5c+)]ra' , 

go = -c+[27c^(c+ - 1) - 21c^c+ + 

4 + c_4(7 + 5c+)]r2a'2 , 

r = 12[4(c+ - 1) + c_(c+ - 2)c+ - 4]m2 X 

{2[-3c_(c+ - 1) + c+]A^a2 + 

[3c-{c+ - 1) + c+]{N\'^ + 1)} . 

The equation for a is 



(A4) 



Y 



(A3) 



where 



gs = -c+[34(c+ - 1) + c_(c+ - 2)c+ + 4]iV x 

[(3c_ - c+)iV2 + 2(3c_ + c+)rNN' + (3c_ - c+yN'^] , 

57 = 2c+[3c?.(c+ - 1) + c_(c+ - 2)c+ + 4]riV2 ^ 

[(3c_ + c+)N + (3c_ - c+)rN']a' , 

ff6 = -{12[3c3 (c+ - 1) + c_(c+ - 1)4 + 4 + clc+ X 

(4c+ - 5)]b27V + [-4(12 + c+) + c_4(12 - 19c+ + 

34) + 9ci(4 - 9c+ + 54) + 3clc+{20 - 33c+ + 84)]^^ 

+4[-34 + c_4(3 - 9c+ - 24) + 9ct{l - 3c+ + 24) + 

3c?_c+(5 - 12c+ + 4:cl)]rNN' + c+[27c^_{c+ - 1) - 

2l4c+ + 4 + c_4(7 + 5c+)]r2iV'2 + c+[-9ci(c+ - 1) + 

3c^c+ + c_(-5 + c+)4 + 4]r2iV3a'2} , 

35 = -2c+rN {[27 c^_{c+ - 1) - 4 - c_4(13 + 5c+) + 

3c^c+(-13 + 6c+)]iV + 2[18ci(c+ - 1) + 3c?_(-5 + 

c+)c+ + 4 + c_4(4 + 3c+)]riV'}a' , 

54 = 12[9ci(c+ - 1) - 4 + c_4(-ll + 3c+) + c^c+ X 

(-19 + I2c+)]B^ + [4(12 + c+) + 9ci(12 - 19c+ + 

74) + c_4(132 - 65c+ + 94) + 3c^c+(76 - 79c+ + 

164)]^^ + 2[-(c+ - 6)4 + c_4(66 - 31c+ - 54) + 

27ci(2 - 3c+ + 4) + 3c2_c+(38 - 37c+ + 64)]^^' - 

c+[45cl(c+ - 1) + 34 + 3c?.c+(-13 + 4c+) + 

c^cl{9 + 7c+)yN^a'^ , 

gs = -2c+r{[94(c+ - 1) + 4 + c-4(l + c+) + 

3c^c+(-3 + 2c+)]A^ + [-27c3 (c+ - 1) + 4 + 

3c?.c+(5 + 2c+) + c_4(-ll + 5c+)]rN'}a' , 

(?2 = -c+{27cl{c+ - 1) + 4 + c_4(-17 + 5c+) + 

34c+(-15 + 8c+) + [-63ci(c+ - 1) + 34 + 

34c+(ll + 4c+) + c_4(-27 + llc+)yNa^} , 

£,1 = -2c+[-27c3_(c+ - 1) + 3c2_(13-6c+)c+ +4 + 



c+r{-l + Na^)Y ' 
where 

/i3 = 2c_4^^[(3c- - c+)N'^ + 2(3c_ + c+)rNN' + 
(3c_ - c+)r2Ar'2] ^ 

/12 - 4c_4riV4p^_ ^ ^ (3^_ _ c+)rAr>' , 

/ii = 4iV2{-24c_c+(c_ + c+)B^N + [9ci(c+ - 
1) - 74 + c_4(17 + 7c+) - 3c2_c+(-5 + 8c+)]N^ + 
2[94(c+ - 1) + c_(13 - 5c+)4 + 74 - 3c^c+(l + 
4c+)]rA^A^' + [9c3_(c+ - 1) + 3^(5 - 12c+)c+ - 74 + 
c_4(17 + 3c+)]r2iV'2 + 2c_(3c_ - c+)clr^ N'^a'^} , 
/lo = 4riV2{[94(c+ - 1) + 154 - 3c2_c+(9 + 4c+) - 
c_4(3 + 5c+)]N^ - 16c+[-2c-c+ + 4 + c^(-3 + 
6c+)]rArAr' + [-9c^(c+ - 1) + 3^0+ + 
c_(-5 + c+)4 + 4]r2iV'2}a' , 

/g = c+7V(-12[3c'l(c+ - 1) + c_(3 - 7c+)4 + 34 - 

cic+{3 + 4:C+)]B'^N + 2[(18 - 13c+)4 + Sclc+{-6 + 

7c+) + c_4(18 - 37c+ + 34) - 9cl(2 - 5c+ + 

34)]?-^^' + 6c+[-3cl(c+ - 1) - 7c_4 + 4 + 

cic+{-5 + 9c+)yN'^ ~ 3c+[3c^(c+ - 1) + 

c_(-5 + c+)4 + 34 + c^c+(-ll + 20c+)]r2iV''^a'2 - 

2N^{2[9ct{c+ - if - 4(3 + 5c+) + c_4(3 + 10c+ + 

24) - 3c^c+(-5 + C+ + 34)] - c+[-9ci(c+ - 1) + 

3clc+ + c_ (-5 + c+)4 + ciyN'a'^}) , 

h = c+riVa'{12[3ci(c+ - 1) + c_(c+ - 1)4 + 4 + 

c^_c+{Ac+ ~ b)]B'^N - 2[4(-6 + 7c+) + c_4(6 - 



23c+ 



9c^(2 - 5c+ + 34) + 3c2_c+(10 - 25c4 



124)]7V2 - 12[3c3_(c+ - if - 4 + c_4(l + 10c+) + 
c2_c+(5 + 4c+ - 134)]7-^^' + 4c_c+[9c2_(c+ - 1) - 
6c_c+ + 4(3 + c+)]r2Ar'2 + c+[-9cl(c+ - 1) + 
ic^_c+ + c_ (-5 + c+)4 + 4]r2iV3a'2} , 
/7 - c+(24[3c3_(c+ - 1) + c_(3 - 4c+)4 + 34 - 
c2_c+(3 + c+)]B2Ar + 2c+[27c3_(c+ - 1) + 4(-48 + 
5c+) + c_c+(-96 + 47c+ + c\) + 3c2_(-16 + 5c+ + 
44)]riVA^' + [4(12 + c+) - c_4(12 - 13c+ + 4) + 
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9c?. (-4 + 3c+ + 4) ~ 3c?.c+(20 - 13c+ + 8c\)yN''^ - 
4c+[-9cl(c+ - 1) + 2c\ + c_4(29 + c+) - 
6c^c+(-6 + 7c+)]r2iV3a'2 + 2N'^{~9c\{2 + c+) + 
c-4(-18 + 21c+ + 4) + 9c^(2 - 5c+ + 34) - 
34c+(-6 + 5c+ + 44) + c+[45ci(c+ - 1) + 4 + 
3c2_c+(-ll + 2c+) + c_4(13 + bc+)]r^ N' a'"^}) , 
fe = c+ra'{-48c_[3c^(c+ - 1) + (c+ - 3)4 + 
2c_c+(-3 + 2c+)]B^N - 4c+[-27c^(c+ - 1) + 
64(10 - 7c+)c+ + 44 + c_4(37 + c+)]iV^ + 
8[4(3 + 2c+) + 9c^(-l + 4) + c_4(-3 + 8c+ 
+24) - 3c^c+(5 - 2c+ + 34)]r7V7V' - c+[27cl{c+ - 
1) - 21c?. c+ +cl+ c_c^(7 + 5c+)]r2A^'2 + 
2c+[27c?.(c+ - 1) + c^ + 3c2_c+(-7 + 2c+) + 
c-cl{7 + 3c+)yN\'^} , 

h = -2(l8[(c+ - 4)c3 + c2_c+(-12 + 7c+) - c_c^ x 
(12 - 5c+ + c^) + ct (-4 + 3c+ + c^)]S2 + 
[-c^(-72 + 30c+ + c\) + 9ci(8 - 10c+ + c\ + c%) + 
c_c^(216 - 150c+ + 23c+ + c\) + 3c^c+(72 - 70c+ + 
llc^ + 2c%)]rN' + c+[(6 - 5c+)c3 + c_c^ (66 - 
119c+ - 7c\) + 27c3 (2 - 3c+ + c\) + 3c^c+(38 - 
65c+ + AQc\)yN'^a''^ + A^{-2c3 (-36 + 9c+ + c\) + 
184(4 - 5c+ + ci^) - c_c^(-216 + 126c+ - 16c^ + 
c\) - 3c^c+(-72 + 66c+ - 6c^ + ci^) + c^[634(c+ - 

I) - 45c^c+ + c\ + c_c^(19 + c+)]r^ N' a'^}) , 

fi = -2c+ra'{-6[9c^(c+ - 1) - c^ + c_c^(3c+ - 

II) + c2_c+(-19 + 12c+)]B2 + [(6 - 7c+)c3 + c_c^ x 
(66 - 73c+ - 7c;^) + 9ci(6 - llc+ + bc\) + 3c^c+ x 
(38 - 55c+ + 22c^)]iV + 6[-c3 + c2_c+(-19 + 8c+) + 
c_c^(-ll + 2c+ + c\) + 3c^(-3 + 2c+ + c\)]rN' + 
2c_c+[27c2_(c+ - 1) - 18c_c+ - (-9 + c+)c\yN^a'^} , 



/g = c+{-36c^ - 84c^c+ + 27c^c+ - 60c_c^ + 
39c^c^ + 9c^c^ - 12c^ + 13c_c^ + c| - c_c| + 
2[4c3 (3 + c+) + 18c3 (6 - 7c+ + c\) + c_c^(132 - 
110c+ + 3c^) + 3c2_c+(76 - 80c+ + ilcj\)yNa''^ - 
2c+[-274(c+ - 1) + c^ + 3c^c+(5 + 2c+) + 
c_c^(-ll + 5c+)]r3iV'a'2} , 
/2 = c+ra'{cX + c_c^(144 - 41c+ - 7c^) + 
9c?.(16 - 19c+ + 3c^) + 3c^c+(96 - 71c+ + 12c^) - 
2c+[-45ci(c+ - 1) + c^ + 3c2_c+(9 + 2c+) + 
c_c^(-17 + 3c+)]r2iVa'2} , 

/i - -c+[(-36 + c+)c\ + 9c3_(12 - 13c+ + c\) + 

c_c^(36- 35c+ + llc^) + 
3c^c+(60 - 51c+ + 2Qc\)]r^a^ , 

/o = -c^[27c^(c+ - 1) - 21c^c+ +c\ + 

c_c^(7 + 5c+)]r3a'3 . 

If we remove N' from (jA3p and (jA4[) . we can write 
them as the form in (|3.6p and p.7p . Since it is too te- 
dious, we do not perform it. From (|A2[) to (jA4[) . we 
obtain the asymptotic form for r ^ 00 as 



N{r) 


= B{oof 


Ni 

+ — + •••, 

r 


(A5) 


B{r) 


= B{oo) -f- 




(A6) 


a{r) : 


1 

" B(oo) ^ 


r 


(A7) 



where iVi, Bi and ai are constants. 

We should be careful about (— 1 + Na^) in the de- 
nominator in (jA4[) since (|3.9p and (|3.1ip show that 
(— 1 + No?') asymptotically approaches zero. However, 
since (—1-1- Na^) oc 1/r for r — s- cx), this is canceled by r 
in the denominator in (|A4p . 
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